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Bj- 

1 Introduction 

■ Differential equations defined on a complex domain frequently appear in mathematics and 
I/-) , physics. One of the most important differential equations is the Gauss hypergeometric differen- 
tial equation. Mathematically, it is a standard form of the second-order differential equation 
with three regular singularities on the Riemann sphere. Global properties of the solutions, i.e., 

f"*- . the monodromy, often play decisive roles in applications to physics and mathematics. 

A canonical form of a Fuchsian equation with four singularities is given by Heun's equation 
which is written as 

> ■ 

d adw — q \ ~, . . 

T+ ( iv *\ /H = °» L1 
aw w{w — l)(w — t) I 




with the condition r y + 5 + e = a + /3+l. This equation appears in several topics in physics, 
i.e., astrophysics, crystalline materials and so on (see [37] and references therein). Although the 
problem of describing the monodromy of Heun's equation is much more difficult than that of the 
hypergeometric equation, Heun's equation has been studied from several viewpoints. A method 
of finite-gap integration is available on a study of Heun's equation, and consequently we have 
some formulae on the monodromy. 

2g-l 

If there exists an odd-order differential operator A = (d/dx) 2d+1 + ^ bj(%) (d/dx) 29 J 

j=0 

such that [A, —d 2 /dx 2 + q(x)] = 0, then q(x) is called the algebro-geometric finite-gap potential. 
Note that the equation [A, —d 2 /dx 2 +q(x)] = is equivalent to the function q(x) being a solution 
to a stationary higher-order KdV equation (see [TO]). In our setting, finite-gap integration is 
a method for analysis of the operator —d 2 /dx 2 + q(x) where q(x) is an algebro-geometric finite- 
gap potential. Originally, the finite-gap property is a notion related to spectra. Let H be the 
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operator —d 2 /dx 2 + q(x), and the set <Jb{H) be defined as follows: 

E G (Jb(H) 44> Every solution to [H — E)f(x) = is bounded on x G 
If the closure of the set (Jb(H) can be written as 



a b (H) = [E , Eh] U [E 2 ,E 3 ] U • • • U [E 2g , oo), 

where i?o < £-1 < • • • < E 2g , i.e., the number of bounded bands is finite, then q(x) is called 
the finite-gap potential. It was established in the 1970s that, under the assumption that q(x) is 
a periodic, smooth, real function, the potential q{x) is finite-gap if and only if q(x) is algebro- 
geometric finite-gap. 

On the approach by the finite-gap integration for Heun's equation, it is essential to transform 
Heun's equation into a form with the elliptic function. The transformed equation is written as 

(H - E)f(x) =(-^ + J2 k{k + l)p(x + Ui) - e] f{x) = 0, (1.2) 



i=0 



where p(x) is the Weierstrass p-function with periods (2ui, 2los), uq(= 0), u>i, uj 2 {= —oj\ — us), 
W3 are half-periods, and li (i = 0,1,2,3) are coupling constants. Here the variables w and x 
in equations (|1.2jl are related by w = (p(x) — p(ui)) / (p(u 2 ) — p{ui{)). For details of the 

transformation, see [34^ [38j 143] . The expression in terms of the elliptic function was already 
discovered in the 19th century, and some results which may relate to finite-gap integration were 
found in that era. For the case when three of Iq, l±, l 2 , I3 are equal to zero, equation (|1 .2j) 
is called Lame's equation. Ince |19] established in 1940 that if n G Z>i, u\ G R \ {0} and 
UJ3 G \J — 1R \ {0}, then the potential of Lame's operator 

d 2 

~dx 2 + U<yU + ^ P<yX + W3 ' ) ' 

is finite-gap. In the late 1980s, Treibich and Verdier [48] found that the method of finite-gap 
integration is applicable for the case lo,h,h,h £ Z. Namely, they showed that the potential in 

equation (II. 2h is an algebro-geometric finite-gap potential if Zj G Z for all i G {0, 1,2,3}. There- 

3 

fore the potential Yl k(k + ^)p{ x + ^i) is called the Treibich- Verdier potential. Subsequently 

i=o 

several others [T71 [3S1 H21 IIH 113 US] have produced more precise statements and concerned 
results on this subject. Namely, integral representations of solutions [381 142j. the Bethe Ansatz 
|17| 142] . the global monodromy in terms of the hyperelliptic integrals [44] , the Hermite-Krichever 
Ansatz [35] and a relationship with the Darboux transformation [36] were studied. 

In this paper, we discuss some approaches to finite-gap integration for multivariable cases. 
A multivariable generalization of Heun's equation is given by the Inozemtsev model, which 
is a generalization of the Calogero-Moser-Sutherland model. The Inozemtsev model of ty- 
pe .BCjv [20j is a quantum mechanical system with A^-particles whose Hamiltonian is given by 



d 2 

j=l 3 l<j<k<N 
N 3 

+ Y + l )P( X 3 + 

j=l i=0 



where / and l{ (i = 0,1,2,3) are coupling constants. It is known that the Inozemtsev model 
of type BCn is completely integrable. More precisely, there exist operators of the form = 
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N 



2k 



£ 3§- + (lower terms) (k = 2, . . . , N) such that [H, H k ] = and [H kl , H k2 ] = (k, fa , k 2 = 
j= i V J / 

2, . . . ,N). Note that the Inozemtsev model of type BCn is a universal completely integrable 
model of quantum mechanics with Bjy symmetry, which follows from the classification due 
to Ochiai, Oshima and Sekiguchi [30]. For the case N = 1, the potential coincides with the 
Treibich-Verdier potential and the spectral problem for the Inozemtsev model of type BC\ is 
equivalent to solving Heun's equation. 

By the trigonometric limit (t(= u^/loi) — > \J — loo), we obtain the trigonometric Calogero- 
Moser-Sutherland model. The trigonometric model is well-studied by multivariable orthogo- 
nal polynomials (i.e., the Jack polynomial and the multivariable Jacobi polynomial). Vadim 
Kuznetsov and his collaborators studied multivariable orthogonal polynomials from the aspects 
of separation of variables [25], the Pfaff lattice [1] and the Q operator [23]. Note that in the 
paper [26] relationships among separation of variables, integral transformations and Lame's 
(Heun's) differential equation were discussed. Applications of these technique for models with 
elliptic potentials are anticipated. 

Although the Inozemtsev model of type BCn is much more difficult than the trigonometric 
one, some approaches (perturbation from the trigonometric limit, quasi-solvability etc.) were 
introduced. Now we hope to develop analysis of this model by finite-gap integration. Although 
we can regard several works to be on this direction, in my opinion, they are still far from complete 
understanding of the model. We may consider the multivariable Darboux transformation as 
a possible approach, but we should develop it in the future. 

This paper is organized as follows. In Section [21 we review the finite-gap integration of Heun's 
equation. An approach by the Darboux transformation is introduced. In Section [3l we collect 
results on the Calogero-Moser-Sutherland model and the Inozemtsev model. In Section U] we 
discuss some approaches to finite-gap integration for those models. 

2 Finite-gap integration of Heun's equation 
2.1 Darboux transformations and Heun's equation 

We consider the finite-gap property of Heun's equation in the elliptic form. It is known that 
the Treibich-Verdier potential is algebro-geometric finite-gap, i.e., there exists a differential 
operator A of odd order which commutes with the operator fl"(Wi>W3) where 



In this subsection, we construct an odd-order differential operator A by composing the Darboux- 
Crum transformation which we will explain below. 

We review the Darboux transformation. Let <f>o(x) be an eigenfunction of the operator H = 
—cP/dx 2 + q(x) corresponding to an eigenvalue Eq, i.e. 



For this case, the potential q(x) is written as q(x) = (4>q(x) /4>o(x))' + (4>'q{x) / '4>q{x)) 2 + Eq. 
By setting L = d/dx — 4>' (x) / 4>q(x) and LA = —d/dx — 4>q(x)/(/)q(x), we have the factorization 
H -E Q = L^L. We set H = W + E . Then we have H = -d 2 /dx 2 + q(x) - 2(0 o (x)/0 o (x)) / 
and the relation HL = LH. Hence, if (j)(x) is an eigenfunction of the operator H corresponding 
to the eigenvalue E, then L<j){x) is an eigenfunction of the operator H corresponding to the 
eigenvalue E. This transformation is called the Darboux transformation. We generalize the 
operator L to be the differential operator of higher order in the following proposition. 




(2.1) 
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Proposition 1 (cf. [3]). Suppose that the operator H = —d 2 /dx 2 + q(x) preserves an n-dimen- 
sional space U of functions. Let L be the monic differential operator of order n which annihilates 
all functions in U , and write 

'-(£)"+X>m(£ 

Set 

H = -d 2 /dx 2 + q(x) + 2c' 1 (x). 
Then we have 
HL = LH. 

We call the operator L in Proposition Q] the generalized Darboux transformation or the 
Darboux— Crum transformation. For the case n = 1, let 4>o( x ) be a non-zero function in U, then 
U = C4>q(x), the operator which annihilates (f>o(x) is given by L = d/dx — 4>' (x) / <f>o(x) and 
the operator H is given by H = H — 2(c/)' Q (x) / <Pq{x))' . Hence the proposition reproduces the 
Darboux transformation for the case n = 1. 

We apply Proposition[T]to Heun's equation. For this purpose, we recall the quasi-solvability of 
Heun's equation. If a finite-dimensional space is invariant under an action of the Hamiltonian H 
of a model, the model is partially solvable by diagonalizing the matrix representation of H. This 
situation is called quasi-solvable. On the quasi-solvability of Heun's equation, we have 

Proposition 2 ([<|3], Proposition 5.1). Let cti be a number such that on = —l{ or ai = Zj + 1 

3 

for all i G {0, 1, 2, 3}, and set d = — ^ «j/2. Suppose that d G Z>o, and let V^ 0iQ1|Q . 2|Q . 3 be the 

i=0 

d+ I -dimensional space spanned by 



{$(p(x))p{x) n } 

k ) n=0,. 



where 8(z) = [z - ei) ai/2 (z - e 2 ) a2,2 {z - e 3 ) Q3/2 . Then the operator ffCWi.Ws) f see eg?ia . 
preserves the space V a0)ai)a2j a 3 - 

By applying Propositions [Hand [21 we obtain the following proposition after some calculations: 

Proposition 3 (|46j). Let on be a number such that on = —li or on = k + 1 for all i G {0, 1, 2, 3} ; 

3 

and set d = — Q i/2- Suppose that d G Z>o, and let L a0tO , ljOl2jOl3 be the monic differential 
operator of order d + 1 which annihilates the space V a0jaii a 2tO , 3 . Then we have 

TT(a +d,ai+d,a 2 +d,a 3 +d) t _ j fj(lo,h,h,h) 

11 -^00,01,02,03 Oo,Ol, 02,03-*-* 

We construct an odd-degree differential operator A which commutes with H by composing 
the operators L ao>ai>ol2t0 , 3 . We define the operator L a0jaijOl2tO , 3 as follows: 

3 

-^Oo, Q l, Q 2,a3 ) G Z<o; 

1=0 

I J 3 

-^00,01,02,03 \ r In <- n, 

-^1-q ,1-oi,1-«2,1-O3, 2^ a i/ z G £i>2, 

i=0 

1, otherwise. 
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Set 



lo = (-h + h + h + h)/2, 
1% = (lo + h-h + h)/2, 



l{ = (k -h + h + h)/2, 
q = (l Q + l 1 + l 2 - l 3 )/2, 



V 



(/o + h + h + k + l)/2, l{ = (l + h-h-h- l)/2, 



l°2 = (k -h+h-h- 1)A l° 3 = (/ -h-h + h- l)/2. 
The following proposition is proved by applying Proposition [3] four times: 
Proposition 4 ( |46j). Assume li G Z>o /or i = 0, 1, 2, 3. If 1$ + l\ + I2 + h is even, we set 

A = ^-i|,Z|+l,if-|-l -Zg-^-Zi.Zo+l ,-/3,^2-|-l- C '-/§ -Zf^l+l^l+l-^-io-^l -/2,-^3) (2-2) 
while if lo + h + h + h is odd, we set 

A = Ljo+i-jg _jg ,-i°L_ h -i ,i 3+ i -i 2 L_io jt o +1 _ t o _ t oLi 0+1 _ h _ h _ h . (2.3) 
We then have that the operator A commutes with H "o,h,h,h) } j e > 

It is known that, if lo,h,l2,h G ^>o, then there exist four invariant spaces of the operator 
H(lo,h,h,h) j w hich we consider in Section 12.31 and the four operators in Proposition H] are related 
to the four spaces. 

Let hi be the rearrangement of such that k$ > k\ > &2 > &3(> 0). Set 



f fe , I0 + I1 + I2 + h 

(k + fci + k 2 - k 3 )/2, h + h + l 2 + h 

ko, lo + h + h + h 

{ (k + k 1 + k 2 + k 3 + l)/2, h + h + h + h 



even, ko + k 3 > k\ + k 2 ; 

even, ko + k 3 < k\ + k 2 ; 

odd, fc > &1 + k 2 + A; 3 + 1; 

odd, k < k\ + k 2 + k 3 + 1. 



Then g G Z>o and the degree of the operator A is 2g + 1. 

For the case /q = 2, Zi = £2 = ^3 — 0, we have g = 2 and the operator A is expressed as 



^2,-1,-1,0-^1,-2,1,0-^0,2,-1,-1-^- 



2,0,0,0 



d | 1 p'(x) | 1 p'(x) 



+ 



p'Or) 1 p'(x) 



(ire 2 p(x) — ei 2 p(x) — e2 J \dx p(x) — e± 2 p(x) — e2 
d_ _ p'{x) 1 p'{x) 1 p'{x) 
dx p{x) -e x 2 p(x) - e 2 2 p(x) - e 3 



1 



p'(*) 



+ 



P'{x) 



+ 



P'{x) 



2\p(x)-e 1 p(x)-e 2 p(x) - e 3 J dx 



dx) ^ \dx J 2 ^ ^ ^ \ dx 



9 5p(x) 



A 92 )Tx 



2.2 Application of finite-gap property 

We investigate Heun's equation in the elliptic form 



(H-E)f(x) = 0, H = -^ + u{x) 



u{x 



i=0 



(2.4) 
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by applying the finite-gap integration, which is based on the commutativity of H (= —d 2 /dx 2 + 
u(x)) and an odd-order differential operator A. 

Since A is a monic differential operator of order 2g + 1, it can be expressed in the form 



A=(-l)'X)(a,-(x)^ + 6 i (x)) H°- 

j=0 V J 



where clq(x) = 1. We have 



= [(-lYA,H] = J2 

3=0 



H 9-3 



i=o 

9 



i=o 

9 



^ ( 2a! j (x)(-H + tt(x)) + (a?(x) + 26j-(ac))— + Oj-(a;)u'(x) + 6j(ar) ) 



^ ( (a?(z) + 2&J(x))^ - 2a;. +1 (x) + 2a^(x)u(x) + 5 i (x)«'(x) + b'^x) ) tf^'. 
i=o 



Hence we obtain 

bj(x) = -a' j (x)/2 + cj, a"'(x) - 4u(x)d' j (x) + Ad' j+1 (x) - 2v!(x)dj(x) = 
for some constants Cj (j = 0, . . . , g). Therefore we have 

Proposition 5. Set ao(x) = 1 and d g+ \{x) = 0. The operator A may be expressed in the form 



A =(-!)< 



9 f 

i=o L 



d I f d 



dx 2 V dx 



.(, ))}//" < • £,-,•//•* 



(2.5) 



/or some functions dj(x) (j = 1, . . . , g) and constants Cj (j = 0, . . . , g), where the functions 5,(x) 
(j = 0, . . . ,g) satisfy 



d'j'(x) — 4it(x)aj(x) + 4aj +1 (x) — 2u (x)dj(x) = 



(2.6) 



We define a function 3(x,E) which plays the important role for the solutions and the mon- 
odromy of Heun's equation. Set 



l(x,E) = J2a g _ i (x)E i . 



(2.7) 



i=0 



It follows from equation (|2.6|) that E(x,E) satisfies a differential equation satisfied by products 
of any pair of the solutions to equation (|2.4D . i.e., 



(2.8) 



{£s-Mn(x)-E)^--2u'(x))E(x,E)=0. 



On the basis of Proposition [5] and the function S(x, E) in equation (|2.7|) . we have 
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Proposition 6 ([46]). (i) The constants cj (j = 1, . . . ,g) in equation (|2.5p are all zero, 
(ii) The function E,(x,E) is even doubly-periodic and expressed as 

3 h-i 

Z(x,E)=c (E) + ^r / ^Tbf\E)p(x + u; i ) l >-i, (2.9) 

i=0 j=0 

where the coefficients cq(E) and bj(E) are polynomials in E. The coefficients do not have com- 
mon divisors and the polynomial cq(E) is monic. We have g = deg E cq(E) and the coefficients 
satisfy deg E b^(E) < g for all i and j. 

For the case lo = 2, h = l 2 = h = 0, we have 

S(x, E)=E 2 + 3Ep(x) + 9p(x) 2 - \g 2 . 

where = p(ui) and g 2 = -4(eie 2 + e 2 e 3 + e 3 ei). 

Note that the function S(x, E) can be also obtained as the function satisfying equation (j2.8|) 
and Proposition [6] (ii) (see [4*2]). 

We can derive an integral formula for a solution to equation (|2.4j) in terms of the doubly 
periodic function H(x, E). Set 



Q(E) = E(x, E) 2 [e-Y, k{k + l)p(x + 




\ i=0 
1 , ,d 2 E(x,E) 1 , 

It is shown by differentiating the right-hand side of equation (|2.10p and applying equation (|2.8p 
that Q(E) is independent of x. Thus Q(E) is a monic polynomial in E of degree 2g + 1, 
which follows from the expression for "E(x,E) given by equation (|2.9p . For the case Iq = 2, 
h = h = h = 0, we have 



Q(E) = (E 2 -3g 2 )H(E-3e l ). 



i=l 



The following proposition on the integral representation of a solution to equation (|2.4p was 
obtained in 



Proposition 7 ([42J. Proposition 3.7). Let r,(x,E) be the doubly periodic function defined in 
equation (|2.7p and Q(E) be the monic polynomial defined in equation (|2.10p . Then the function 



A(x,E) = v W^)ewj^^ 

is a solution to the differential equation (|2.4|) . 

If the value E satisfies Q{E) 7^ 0, then the functions A(x, E) and A(— x, E) form a basis of 
solutions to equation (|2.4|) . Since equation (|2.4|) is doubly-periodic, the functions h(x -\-2u)f~,E) 
and A(— (x + 2ujk),E) are also solutions to equation (|2.4p . We consider the monodromy on the 
functions A(x, .E) and A(— x, E). Note that, if A(x + 2u>k, E) is expressed as B(E)A(x, E), then 
A(— (x + 2cok),E) is expressed as B(£') _1 A(— x, E). We will express B(E) as a hyperelliptic 
integral of second kind. We rewrite the function H(x,i?) and define a(E) as follows: 



l(x,E) =c(E) + J2T,a?(E) (f) 2j p(x + ^), a(E) = £ a® (E) . 



(dx) 

i=0 j=0 v 7 i=0 
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Proposition 8 ( [44j . Theorem 3.7). Assume li € Z>o (i = 0, 1,2,3). Let Eq be a value such 
that Q(E ) = 0. Then A(x + 2cv k ,E ) = (-l) 9fe A(x, S ) /or % G {0, 1} (k = 1, 3) and we have 



A(x + 2w fc , £) = (-l) 9fc A(x, E) exp 



1 



-2rjka(E) + 2w fc c(£) 



where g k = C( w fc) = 1, 3) and C( x ) ^ s the Weierstrass zeta function. 



For the case Iq = 2, l± = l 2 = h = 0, we set 
function a(E) and c(i£) are determined as 



/ Sg2- Then (?i = 93 = and the 



c(E) =E 2 - \g 2 , 



Hence we have 



a (E) = 3E. 
( 



A(x + 2u k , E) = A(x, E) exp 



V 



E -6E Vk + (2E 2 - 3g 2 )iv k 
J-(E 2 -3g 2 )ll(E-3e i 



-AE 



I 



We review the propositions related with the Bethe Ansatz (Proposition [9]) and the Hermite- 
Krichever Ansatz (Proposition I lOj) . which are also reductions of the finite-gap property. 

Proposition 9 ([42] . Theorem 3.12). (i) If the value E satisfies Q(E) 7^ 0, then there exists 
t\, . . . ,t n and C such that tj 7^ tji (j 7^ j'), tj <jj{L + uj^TL and A(x, E) is expressed as 



A(x,E) = C 



11 (r(x + tj) 



0-(xy°0-l(x) h O-2{x) l2 0-3(x) l 3 



exp — x 



(2.11) 



where o~(x) is the Weierstrass sigma function and o-j(x) {i = 1,2,3) are the co-sigma functions 
defined by 

Ui(z) = exp(-r]iz)a(z +Ui)/a(uji). 

(ii) The function 



A(x) 



a(x) lo ai(x) ll a2{x) l2 a3(x) 13 



exp(cx), 



(2.12) 



with the condition tj 7^ ty (j 7^ j') and tj u)\L + uj-^L is an eigenfunction of the operator H 
(see equation (|2,4p ). if and only iftj (j = I, . . . ,1) and c satisfy the relations, 

3 

2 + **) - foCHi) - £ h(({-tj + Wi) - C(w0) = -c, (j = 1, . . . , /), (2.13) 



1=1 



(l-<S, ,o) c + ^Cfe) =0, 



1 - *,. i0 ) [ c + lC(ui) + E C(-Wi + tj) ]=0, (i = 1, 2, 3). 
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The eigenvalue E is given by 

3 

E = -c 2 + (l h + hh)ei + (l h + hh)e 2 + {kh + hh)e* ~ Yl l ^ 2c + ^) 

i=l 

I 3 

j=l i=0 j<k 

Equation (|2.13|) is called the Bethe Ansatz equation for the Inozemtsev model of type BC\ 
(see [12]). Note that Gesztesy and Weikard [T7] obtained similar results. The monodromy of 
the function A(x) in equation f|2. 12j) is written as 

A(x + 2co k ) = exp(2? ?fc (i 1 + • • • + t t ) + 2co k (c - ((h) C(*l)))A(z) 

for k = 1,2,3. 

In order to describe the proposition on the Hermite-Krichever Ansatz, we define 

T / \ a(x + u>i — a) />./ \ \ 
$i x, a) = y 1 exp C a )x , i = 0, 1, 2, 3 . 

cr(x + CJi) 

Proposition 10 ([45J). There exist polynomials P\(E), . . . , Pq(E) such that, if P2(E) ^ 0, then 
A(x, E) is written as 



A(x, £7) = exp (kx) £ £ ( — <^(x, a) (2.14) 



for some values b$ (i = 0, . . . , 3, j = 0, . . . , Zj — 1), a and k. T/ie values a and k are expressed 
as 

For i/ie periodicity of the function A(x,E), we have 

A(x + 2wfc, £7) = exp(-2?? fc a + 2u) k C,(a) + 2«a;fc)A(x, E) 
for k = 1, 3. 

Note that q in equation (|2.14p and tj in equation (|2.1ip satisfy the relation a = — tj- To 

3=1 

calculate the polynomials P\ (E),... ,Pq(E), it is effective to apply the notions "twisted Heun 
polynomial" and "theta- twisted Heun polynomial" (see [45]). 

If Iq = 2, l± = 1% = I3 = 0, then the values a and k are expressed as 

(E-3 ei )(E + 6 ei ) 2 2 r —— 

P{a) = 61 W 2 ~ 3g 2 ) ' = 3(^-3fl,r ^ 

2.3 Relationship among commuting operators 

We review a relationship among the operators H, A, the polynomial Q(E) and the invariant 
subspaces. On the operators H and A, we have the following relation: 
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Proposition 11 ( |44j . Proposition 3.2). Let H and A be the operators defined by equation (j2.4j) 
and equations (|2.2|) , (|2,3|) . and Q(E) be the polynomial defined in equation (|2.10|) . Then 

A 2 + Q(H) = 0. 

It is known that, if lo,h,h,h £ ^>0i then there exist four invariant subspaces with respect 
to the action of the operator H. We describe the spaces more precisely. Let V Q0iQliQ2ia3 be the 
space defined in Proposition [2] and 

3 
i=0 

{0}, otherwise. 
If lo,h,hih G ^>o an d /o + ^i + ^2 + ^3 is even, then the operator H preserves the space 

v = u^i ^ h -i 2 -i 3 © U-i -i u i 2+1j i 3+1 © £/_ i0)il+ i ,-{2,13+1 © c/_i 0)Il +i jj2 +i _j 3 , (2.15) 

and also preserves the components in equation (|2.15p . If lo,h,l2,h G ^>o and /o + h + ^2 + ^3 
is odd, then the operator H preserves the space 

V = U-ig-h-hjs+i © U-i -i lt i 2+ i-i 3 © U-i 0t i 1+ i-i 2 -i 3 © Ui +i ,-i 2 ,-i 3 , (2.16) 

and also preserves the components in equation (|2.16p . Then we have 

Proposition 12 ([M|). (i) The operator A annihilates any elements in the space V . 

(ii) The monic characteristic polynomial of the space V with respect to the action of H 
coincides with Q(E). 

Note that the operator A was constructed by composing the generalized Darboux transfor- 
mations which are related to the spaces in the components of equation (|2.15p or equation (|2.16p . 

3 Results on the Calogero— Moser— Sutherland model 
and the Inozemtsev model 

We are going to consider multidimensional generalizations of Lame's equation and Heun's equa- 
tion in the elliptic form. For this purpose, we introduce the quantum mechanical systems. 

3.1 The elliptic Calogero— Moser— Sutherland model 

The elliptic Calogero-Moser-Sutherland model (or the elliptic Olshanetsky-Perelomov mo- 
del [31]) of type An_i is a quantum many-body system whose Hamiltonian is given as follows: 

1 N d 2 

i=l ' l<i<j<N 

where p{x) is the Weierstrass elliptic function. For the case N = 2, the model reproduces Lame's 
equation by setting x\ — X2 = x and restricting to the line x\ + X2 = 0. 

This model is known to be completely integrable, i.e., there exist 7V-algebraically independent 
commuting operators P}. (k = 1, . . . , N) which commute with the Hamiltonian H. Namely, by 
setting 

Pk = E 2iY E a ^ - x ^ x * - x 4 ) • ■ ■ 

0<j<[fc/2] J ' V Jh creSiv 
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x p{x 2 j-i - x 2j )d 2j+1 d 2j+2 ■■■d k ), (3.1) 

where Sn is the symmetric group, [x] is the integral part of x and di = d/dxi, we have H] = 
(1 < k < N) and [Pk,Pk>] = (1 < k, k' < N) (see [30]). The Hamiltonian H is expressed as 
H = P 2 - P?/2. 

By the trigonometric limit (r — ► y 7 — loo) of the elliptic Calogero-Moser-Sutherland model 
where (1, r) is the basic periods of the elliptic function, we obtain (up to an additive scalar) the 
Hamiltonian of the trigonometric Calogero-Moser-Sutherland model, 

1^ 3 2 ^ 1 



= " 2 E ^ + ^/(i + 1) E sin 2 (7r( ))' 

1=1 * Ki<Ti<V V V « J 1 1 



i=l 1 l<i<j<N 

The eigenstates of the Calogero-Sutherland model are described by the Jack polynomial 
J^\x) (A G Mjv) (see [39]), i.e., 

i? trig (jl^ T) (X)A(X) /+1 ) = (e + 27T 2 E[ T ^ ] )/ x T ^ ) (X)A(X) l+ \ 
where Xi = exp (27T\/^Txi) , Mn = {A = (Ai, A2, • • • , Aj\r)|« > J =^ Aj — Aj G Z>o}, eo = 
Ivr 2 (/ + l) 2 iV(iV2-l), A(X) = (X 1 X 2 ...X N )— Ui^-Xj) and sf 1 = £ A 2 +£ ^^A,. 

t<j i=l i=l 

In particular, the ground-state is given by A(X)' +1 . Several properties of the Jack polynomial 
were studied. Vadim Kuznetsov and his collaborators studied the Jack polynomial and related 
polynomials from the aspects of separation of variable [25], the Pfaff lattice pQ and the Q 
operator |24j . 

In contrast with the trigonometric models, the elliptic models are less investigated and the 
spectra or the eigenfunctions are not sufficiently analyzed. There is, however, some important 
progress due to Felder and Varchenko. They introduced the Bethe Ansatz method for the N- 
particle elliptic Calogero-Moser model with the coupling constant I a positive integer. Note 
that Hermite essentially introduced the Bethe Ansatz method for the case N = 2 and I E Z 
(see [M]), and Dittrich and Inozemtsev [9] did it for the case N = 3 and I = 1 in a different 
representation. 

Fix the parameters N and I. We set m = IN(N— 1)/2. Let c : {1, . . . , m} — > {1, . . . , N} be the 
unique non-decreasing function such that c _1 (j) has (N — j)l elements. Let (1 < i < N) be an 

N N 

orthonormal basis of M N with an inner product (•,•). Set ai = e^— ej+i, fj* = Xi€i\ x i = 0}j 

i=l i=l 

Pi = i(2N - i - 1)1/2 and Vi = + \,Vi-\ + 2, . . . , Pi } (1 < i < N - 1). Let W be the 

set of maps w = (w\, . . . , wn) {wi : Vi — > {i, i + 1, . . . , N — 1}) such that #{^ rl (i)} = I for 
1 < i < j < — 1. For w = (wi, . . . , wn^i) £ W, let F w be the set of maps / = (/1, . . . , fN-2) 
(fi : — > V^) such that (i) fa is injective (ii) If /j(x) = y then w i+ i(x) = Wi(y). Set 

00 

\ I 1 \n— 1 / I i"/„ 1 /o\2"l 



6»i(x) = 2 2^(-l) exp (rvrV-Un - 1/2)^) sin(2n - l)vrx, 

n=l 

0i(x) fl'WflQK - A) 

For £ G f)*, we introduce the functions $v(ti, . . . , t m ) and o>(t; x) as follows 
^ T (t 1 ,...,t m ) = e 27T ^ T ^ t ^ ) 

n n e ^-h? n 



X 



l<J<(iV-l)/ =(»)=c(j) |c(i)-c(j)|=l 



12 



K. Takemura 



N-l Pi 

^^e^MJ] EH II ^ m (t*-t A w). 

ujsw/sF™ i=i fc=pi_i+i 
where to = 0, fo(k) = 0. Then we have 

Proposition 13 f [12[ \13\ lllj). // (t?) • • • j tm) satisfy the following Bethe Ansatz equations, 
i/ie function Lo(t°;x) is an eigenfunction of the Hamiltonian H with the eigenvalue 



2vr 2 (£, - 2 7 rv /Z T| : 5(tS, . . . , r) - Z(Z + 1)(N - 1)N V , 



where 



S(h, . . . ,t m ;r) = ^2(a c{i) ,a c{j) ) log B(ti - tj) - ^ lN\ogO(ti 

i<j c(i)=l 



»7 




~ 4 St3^ and P = exp(27r^r). 

n=l ^ / 

Therefore, if we find solutions to the Bethe Ansatz equations, we can investigate the Calogero- 
Moser-Sutherland model in more detail. There are two things to be considered for applying 
Proposition [13] to the spectral problem of the elliptic Calogero-Moser-Sutherland model. The 
first one is to find the condition when the eigenfunctions obtained by the Bethe Ansatz method 
are connected to square-integrable eigenstates and the second one is how the solutions of the 
Bethe Ansatz equation behave. 

On the first question, the condition is described as the parameter £ belonging to some lattice 
(the weight lattice of type Ajy^i). By symmetrizing or anti-symmetrizing the function ui(t°; x), 
we obtain square-integrable eigenstates, although we must check that they are identically zero 
or not. 

On the second question, we consider the solution at p = exp(2ir\/— It) = (the case of 
the trigonometric limit r — ► \J— loo) and look into the behavior where p is near 0, because it 
is hopeful to solve the Bethe Ansatz equations for the trigonometric case in contrast to being 
hopeless directly for the elliptic case. A key tool to connect the trigonomertic solutions to 
the elliptic solutions is the implicit function theorem. Thus we construct the square-integrable 
eigenstates and obtain the main result in [50], which gives a sufficient condition for regular 
convergence of the perturbation expansion. In particular, for the case N = 2, I G Z>i and the 
case N = 3, I = 1, we have convergence of the perturbation series for all eigenstates related to 
the Jack polynomial. 

Note that this idea can be interpreted to consider the elliptic Calogero-Moser-Sutherland 
model by perturbation from the trigonometric Calogero-Moser-Sutherland model. Convergence 
for the general cases was proved in [23] by another method. Namely, by applying Kato-Rellich 
theory, we have convergence of the perturbation series in p for I > and arbitrary N. The 
eigenvalues and the eigenfunctions are calculated as power series by a standard algorithm of 
perturbation. Remark that Fernandez, Garcia and Perelomov [TJJ derived a fully explicit formula 
for second order in p, and Langmann [27[ [28] obtained another algorithm for constructing the 
eigenfunctions and the eigenvalues as formal power series of p, which also gives a formula for all 
orders in p. 

On the Bethe Ansatz for the elliptic Calogero-Moser-Sutherland model, there are some 
problems to be solved. For example, it has not been shown at the moment of writing that the 
eigenfunction u(t , x) written in the form of the Bethe Ansatz is also an eigenfunction of the 
higher commuting operators P3, . . . , Pjy (see also [35]). 
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3.2 The Inozemtsev model 

We now introduce a quantum mechanical system that is a multidimensional generalization of 
Heun's equation in the elliptic form. 

The Inozemtsev model of type BCn [2UJ is a quantum mechanical system with iV-particles 
whose Hamiltonian is given by 

H = -Hq^, +21 ^ 1 + 1 ) Yl (P(xj - x k ) + p{ Xj + x k )) (3.2) 
N 3 
j=l i=0 

where / and li [i = 0, 1, 2, 3) are coupling constants. This is also a generalization of the elliptic 
Calogero-Moser-Sutherland model of type BCn- 

The Inozemtsev model of type BCn is completely integrable, i.e., there exist TV algebraically 
independent mutually commuting differential operators P k (k = 1, . . . ,N) (higher commuting 
Hamiltonians) which commute with the Hamiltonian of the model, and Oshima [32] described the 
commuting operators explicitly. Note that the Inozemtsev model of type BCn (resp. the elliptic 
Calogero-Moser-Sutherland model of type A n) is a universal completely integrable model of 
quantum mechanics with the symmetry of the Weyl group of type Bn (resp. type An), which 
follows from the classification due to Ochiai, Oshima and Sekiguchi [30|, I33J . For the case N = 1, 
the operator (|3.2h appears in the elliptic form of Heun's equation (|2.4p . Therefore the Inozemtsev 
model of type BCn is regarded as a multidimensional generalization of Heun's equation. 

On the trigonometric limit r — ► \J — loo, we obtain the trigonometric Calogero-Moser- 
Sutherland model of type BCn, and we can investigate the Inozemtsev model of type BCn 
by perturbation from the trigonometric model \23\ Wf\ . 

A method of quasi-solvability is available on the Inozemtsev model of type BCn- Finkel 
et al. studied quasi-solvable models in [151 (EL an d they found several quasi-exactly solvable 
many-body systems including the Inozemtsev model of type BCn- We now describe the finite- 
dimensional spaces which are related to the quasi-solvability. The quasi-solvability with respect 
to the Hamiltonian H was established in [16] and reformulated in |41j . 

Proposition 14 ( |16l 141] ). Let a, b% (i = 0, 1, 2, 3) be the numbers which satisfy a G {— I, I + 1} 
and bi G {-k/2,(k + l)/2} (t = 0,1,2,3). Set 



N 3 

$(z)= H (zj-ztfUYlte-ei) 

l<j<k<N j=li=l 



Assume that d = —{{N — l)a + bo + b\ + 62 + 63) is a non-negative integer. Let W^ m be the 
space spanned by 

^(p(x 1 ), P (x 2 ),...,p(x n )) p(xi) m °mp(x 2 r«M---p(x N ) m °w 

such that rrii G {0, 1, . . . , d} for all i. Then we have 
H ■ Wf m C W° ym - 



The quasi-solvability was extended to the commuting differential operators. 
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Proposition 15 ([H], Theorem 3.3). Assume that d = —((N — l)a + &o + b\ + b<i + 63) is a 

non-negative integer. Then P^ ■ W^ m C W^ 7 ™ for k = 1,2, ... ,N, where VF^ ym is the finite- 
dimensional space defined in Proposition Q3] and Pf. are the commuting differential operators 
which ensure the complete integrability. 

By the quasi-solvability, finitely-many eigenvalues and eigenfunctions are calculated by diag- 
onalizing the commuting matrices simultaneously for the case that the assumption of Proposi- 
tion[T3]is true. The eigenfunctions obtained by the quasi-solvability may not be square-integrablc 
in general, although the eigenfunctions for the case that the parameters a, bo, b\ in PropositionfHI 
are chosen as a = I + 1, 60 = (^0 + l)/2 and b\ = (l\ + l)/2 are square-integrable. 

It seems that an explicit expression of the Bethe Ansatz as Proposition [13] for the Inozemtsev 
model of type BCm and corresponding conformal field theory are not known in the moment of 
writing, although Chalykh, Etingof and Oblomkov [6] gave a general recipe for calculating the 
Bloch eigenfunctions. They showed that these are parametrized by a certain algebraic variety 
(the Hermite-Bloch variety) which can be computed. This would lead to a version of the Bethe 
Ansatz for the models including the Inozemtsev model of type BCn, though these Bethe Ansatz 
equations would be rather complicated. For a special case of the BC2 case, this scheme is worked 
out explicitly in (6J § 6]. 

We hope to investigate the Bethe Ansatz for the Inozemtsev model of type BCn to study 
the model in more detail. 

4 Towards finite-gap integration of the Inozemtsev model 

In Section [21 we reviewed the finite-gap integration of Heun's equation, and observed that the 
existence of commuting operator of odd order plays important roles. 

For a multidimensional generalization of finite-gap integration, Chalykh and Veselov in- 

N 

troduced the notion "algebraic integrability". The Schrodinger operator L = — ^ d 2 /dxf + 

i=i 

u(x±, . . . , xn) is called completely integrable [8], if there exist N commuting operators L\ = 
L, L2, . . . , Ln with algebraically independent constant highest symbols (= £? + ••• + Cat)) 

S2 (£),..., sjv(£) (£? = \/—ld/dxj). For example, the Calogero-Moser-Sutherland models are 
completely integrable. On the model of type An, the highest symbols of Pk (see equation (13.111 ) 

are written as ((— \J — l) k /k\) Yl " " The Inozemtsev model of type BCn is also 

ii<ia<— <ifc 

completely integrable. The operator L is called algebraically integrable in the sense of [8], if L is 
completely integrable and there exists one more operator Lq commuting with Lj (i = 1, . . . , N) 
and the highest symbol so(£) °f -^0 takes the distinct values at the roots of the equations 
= Ei (i = 1, . . . , N) for almost all E^. 

On Heun's equation in the elliptic form, if Iq, l\, I2, h are integers, then it is algebraically 
integrable, because there exists a commuting operator A of odd order. Thus we expect that 
algebraically integrable Schrodinger operators also have rich properties. 

Chalykh and Veselov conjectured [7] that the Calogero-Moser-Sutherland model with integral 
coupling constants are algebraically integrable. For the case of type A j\r, Braverman, Etingof 
and Gaitsgory [3] obtained algebraic integrability. More precisely, they established that, if I is 
a positive integer, then the operator 




is algebraically integrable by applying the Bethe Ansatz due to Felder-Varchenko (see Sec- 
tion 13. 1[) and the differential Galois theory. In [6], Chalykh, Etingof and Oblomkov proved 
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that the Chalykh-Veselov conjecture is true and the Inozemtsev model of type BCn (see equa- 
tion (|3.2|) ) is also algebraically integrable, if I, lo, l±, l 2 , h are a h integers. Their method relies on 
results on the differential Galois theory obtained in [1] and the local triviality of the monodromy. 
On an application of the algebraic integrability, the eigenfunctions of the Baker-Akhiezer (Bloch) 
type are considered. We expect further studies for applications of the algebraic integrability on 
the Calogero-Moser-Sutherland models and the Inozemtsev models. 

The explicit expressions of the extra commuting operators were obtained and investigated by 
Oblomkov, Khodarinova and Prikhodsky |29[ l2"Tl [22] for the case I = 1 on the Calogero-Moser- 
Sutherland model of type A 3 and the case I = lo = 1, l\ = l 2 = l 3 = on the Inozemtsev model 
of type BC 2 . On the Calogero-Moser-Sutherland model of type A 3 with 1 = 1, the Hamiltonian 
and commuting operators which guarantee the complete integrability are given as 

H = ~(dl + d\ + 8l)/2 + 2(pi2 + P23 + P3l), 

Pi = di + d 2 + d 3 , 

P3 = did 2 d 3 + 2p 12 <9 3 + 2p 23 di + 2p 3i d 2 , 

(see equation (|3.ip ) where we have used the notations di = d/dxi and py = p(xj — Xj). The 
additional commuting operators are written as 

I12 = (di - d 3 ) 2 (d 2 - d 3 f - 8p 2 3(a! - a 3 ) 2 - 8 P i 3 (c> 2 - d 3 f 

+ 4(pi2 - P13 - Piz){di - d 3 ){d 2 - d 3 ) - 2(p' 12 + p' 13 + 6p 23 )(ai - <9 3 ) 

- 2 (-Pl2 + 6p'l 3 + P23)(^2 - d 3 ) - 2p'/ 2 - 6pi 3 - 6p 23 + 4(p? 2 + p? 3 + P23) 
+ 8(pl2pl3 + P12P23 + 7pi3p23), 

^23 and J31, which are written by permuting the indices. Then any non-symmetric linear com- 
bination of them, e.g., L4 = I\ 2 + 2I 23 would fit into the definition of algebraic integrability 
(see |21j). Explicit expressions of the extra commuting operators for the models which have 
symmetry of the deformed root system A 3 (m) or B 2 (l,m) were also obtained. 

Another possible method for constructing extra commuting operators is the multidimensional 
Darboux transformation, because the commuting operator for the case of Heun's equation is con- 
structed by composing the (generalized) Darboux transformations. Multidimensional Darboux 
transformations were studied from several viewpoints [21 [18, 36, 5j. In [2], based on the existence 
of an explicit eigenfunction of the Hamiltonian H{= H^) with a certain eigenvalue, an alternate 
Hamiltonian H{= H^ N '), matrix valued operators (i = 1, . . . , N — 1) and supersymmetry 
operators Qj + ± j and Q^+i (j = 0, . . . , N — 1) which connect H^' and H^ +1 ' were introduced 
and studied. On the other hand, we know an explicit eigenfunction of the Inozemtsev model of 
type BC N , if the value d(= -((N- l)a + b + h + b 2 + b 3 )) (a G {-1,1 + 1}, E {~k/2, (k + l)/2} 
(i = 0, 1, 2, 3)) in the assumption of Proposition [T3] is equal to zero. Then the alternate Hamil- 
tonian H with respect to H in equation (|3.2p would be written as 

N d 2 

3=1 J j<k 
N 3 

j=l i=0 

In the moment of writing, we do not know an operator L which directly intertwines the oper- 
ators H and H as HL = LH. We expect applications of the multidimensional Darboux trans- 
formation for the analysis of the elliptic Calogero-Moser-Sutherland model or the Inozemtsev 
model. 
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